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ABSTRACT 

The thermal inertia values of Saturn's main rings (the A, B, and C rings and the 
Cassini division) are derived by applying our thermal model to azimuthally scanned 
spectra taken by the Cassini Composite Infrared Spectrometer (CIRS). Model fits 
show the thermal inertia of ring particles to be 16, 13, 20, and 11 Jm~ 2 K _1 s~ 1/2 
for the A, B, and C rings, and the Cassini division, respectively. However, there 
are systematic deviations between modeled and observed temperatures in Saturn's 
shadow depending on solar phase angle, and these deviations indicate that the ap- 
parent thermal inertia increases with solar phase angle. This dependence is likely 
to be explained if large slowly spinning particles have lower thermal inertia values 
than those for small fast spinning particles because the thermal emission of slow 
rotators is relatively stronger than that of fast rotators at low phase and vise versa. 
Additional parameter fits, which assume that slow and fast rotators have different 
thermal inertia values, show the derived thermal inertia values of slow (fast) rota- 
tors to be 8 (77), 8 (27), 9 (34), 5 (55) Jm^K-V 1 ' 2 for the A, B, and C rings, 
and the Cassini division, respectively. The values for fast rotators are still much 
smaller than those for solid ice with no porosity. Thus, fast rotators are likely to 
have surface regolith layers, but these may not be as fluffy as those for slow rota- 
tors, probably because the capability of holding regolith particles is limited for fast 
rotators due to the strong centrifugal force on surfaces of fast rotators. Other addi- 
tional parameter fits, in which radii of fast rotators are varied, indicate that particles 
less than ~ 1 cm should not occupy more than a half of the cross section for the A, 
B, and C rings. 

Key words: Saturn, Rings; Infrared observations; Radiative transfer, Regoliths 



1. Introduction 

This is our third paper analyzing thermal infrared data of Saturn's rings. In Morishima et 
al. (2009; hereafter Paper I), we introduced our new multilayer model, and applied it to ground- 
based data obtained before the Cassini mission. In Morishima et al. (2010; hereafter Paper II), 
we improved our code, and estimated the values of the albedo and the ratio of fast and slow 
rotators using thermal data radially scanned by the Cassini Composite Infrared Spectrometer 
(CIRS). In the present paper, we estimate the thermal inertia values of Saturn's rings using 
thermal data azimuthally scanned by Cassini CIRS including those in Saturn's shadow. 

Saturn's main rings (the A, B, and C rings and the Cassini division) consist of a numerous 
number of icy particles. The range of the particle size deduced from radio and stellar light 
occultations is roughly 1 cm to 10 m (Marouf et al., 1983; Zebker et al., 1985; French and 
Nicholson, 2000; Cuzzi et al., 2009). The ring spectra in ultra-violet (Bradley et al., 2010), 
near-infrared (Poulet et al., 2003; Nicholson et al, 2008), and far-infrared (Spilker et al., 2005) 
light indicate that surfaces of ring particles are covered by tiny regolith particles ranging from 
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micron to cm sizesQThis means that the largest regolith particles are comparable to the smallest 
free floating ring particles. The abundance of free floating micron-size grains is likely to be 
negligible (Dones et al., 1993), except when spokes form (e.g., D'Aversa et al., 2010). Most 
likely, most of small grains (< cm) stick to larger particles due to the surface energy (e.g., Choksi 
et al., 1993). The composition of ring (or regolith) particles is mostly crystalline water ice 
and the mass fraction of contaminants (Tholins, PAHs, or nanohematite) is less than a percent 
(Poulet et al., 2003; Cuzzi et al., 2009). 

Besides the spectroscopy, thermal infrared observations can provide unique information 
about the surface layers of icy particles and about the particle size. The ring temperature drops 
during particles pass through Saturn's shadow and increases after they exit from the shadow to 
the sunshine. The degree of the temperature drop increases with decreasing thermal inertia of 
the particle, if the particle size is larger than the thermal skin depth for eclipse. The thermal 
inertia, which is given by the square root of the product of the physical density, the specific 
heat, and the thermal conductivity, represents how much the particle suppresses its temperature 
change against illumination flux change. While the uncertainties in the density and the specific 
heat may be only by a factor of 2-3, the thermal conductivity varies over several orders of mag- 
nitude depending on structures of surface layers. A numerical model of Shoshany et al. (2002) 
showed that the thermal conductivity decreases by a factor of 10 4 with increasing porosity from 
to ~ 80%. There has been no direct estimation of the porosity of ring particles. However, the 
internal density of ring particles of ~ 450 kg m -3 or less is favorable in dynamical simulations 
for formation of wakes in the A and B rings (Salo et al., 2001; Stuart et al., 2010), because 
signatures of wakes become too strong in simulations using the density of the solid ice with 
no porosity (~ 900 kg m™ 3 ). This indicates that the porosity of ring particles is likely to be 
larger than 50% and the surface layers' porosity may be even higher. An importance of a net- 
work of microcracks was also pointed out by Kouchi et al. (1992), who found very low values 
of the thermal conductivity of amorphous and crystalline ices (thermal inertia values less than 
10 Jm _2 K _1 s -1 '' 2 ) in their laboratory experiments even though the porosity was low. Overall, 
the thermal inertia value is a good indication for the structure of regolith layers represented by 
poorly known quantities such as the porosity. If the particle size is smaller than the thermal skin 
depth, the interior of the particle is roughly entirely isothermal. In this case, the temperature 
drop in the shadow increases with decreasing particle size and is rather independent of the ther- 
mal inertia. If the particle size is comparable to the thermal skin depth, the temperature drop 
depends on both thermal inertia and particle size. 

Prior to the Cassini mission, the thermal inertia values were estimated using thermal data 
mostly from ground-based observations. Morrison (1974) and Froidevaux et al. (1981) esti- 
mated the thermal inertia of the B ring particles to be > 40 Jnr 2 K _1 s~ 1,/2 and ~ 13 Jm~ 2 K _1 s~ 1/2 , 
respectively, using their ground-based data and applying a simple cooling and heating model 



The estimated size of regolith particles increases with wavelength. It may indicate the regolith particle size 
increases with depth in a single ring particle, as longer waves reach to deeper in a ring particle. On the other hand, 
the estimated size may be biased by observed wavelengths even though regolith particles with different sizes are 
well mixed regardless of depth. 
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developed by Aumann and Kieffer (1973). They also found that particles less than about a cen- 
timeter should not occupy a significant fraction of the ring cross section. Ferrari et al. (2005) 
applied their thermal model to their ground based data and estimated the thermal inertias of the 
B and C rings to be 5+{* JnT 2 K _1 s~ 1/2 and 6!' 2 Jnr 2 K _1 s~ 1/2 , respectively. Leyrat et al. (2008a) 
analyzed both their ground-based observation data and Voyager IRIS data for the C ring, using 
their new monolayer model (Ferrari and Leyrat, 2006), and confirmed the very low thermal in- 
ertia estimated by Ferrari et al. (2005). The thermal inertia values for the A ring and the Cassini 
division have not been estimated so far using ring temperatures in Saturn's shadow. 

Since the Saturn orbit insertion of the Cassini spacecraft in July 2004, Cassini CIRS has 
obtained more than a million spectra (7 jum - 1mm) of Saturn's rings (Flasar et al., 2005; 
Spilker et al., 2005, 2006; Altobelli et al, 2007, 2008, 2009; Leyrat et al., 2008b; Flandes et al., 
2010). In Paper II, we estimated the albedo and the fraction of fast rotators with a good radial 
resolution, applying our multilayer model developed in Paper I to the radially scanned CIRS 
data selected in Spilker et al. (2006). Although the thermal inertia values were also estimated in 
Paper II, the data sets were not very suitable for accurate measurements of the thermal inertia. 
For optically thick rings, the thermal inertia was constrained from the temperature difference 
between the lit and unlit faces, but there was degeneracy between the thermal inertia and the 
fraction of particles bouncing at the midplane. For optically thin rings, the thermal inertia was 
constrained from the temperature difference between the morning and evening temperatures, 
but the accuracy of the estimated thermal inertia was not good because the azimuthal locations 
of the data points were far from Saturn's shadow. In the present paper, we analyze azimuthally 
scanned CIRS data including those in Saturn's shadow and precisely estimate the thermal inertia 
values for Saturn's main rings (the A, B, and C rings and the Cassini division), extensively using 
our multilayer model. The azimuthally scanned CIRS data have advantages in their spatial and 
spectral resolutions, azimuthal coverages both inside and outside Saturn's shadow, and solar 
phase coverages as compared with previous ground-based and spacecraft (Pioneer and Voyager) 
observations. Some of azimuthally scanned CIRS data obtained in the early Cassini mission 
have been already published by Leyrat et al. (2008b). In addition to these data, we use some 
more data obtained later. 

In Sec. 2, we discuss how the thermal relaxation time is related to the thermal inertia and 
the particle size. In Sec. 3, we explain data selections and fitting procedures and show fitted 
parameter values. In Sec. 4, we discuss implication from estimated thermal inertia values and 
future works. A summary is given in Sec. 5. 



2. Thermal relaxation time and thermal inertia 

The physics of eclipse cooling is summarized as follows. When the illumination on a ring 
particle changes with a typical frequency a>, its temperature is modified in a surface layer of 
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characteristic thickness (i.e., thermal skin depth) given as 



, ^ CD 

where K, p, and C are the thermal conductivity, the physical density, and the specific heat of 
the particle, respectively. For the case of eclipse cooling of Saturn's rings, the time scale of 
the cycle, 2n/co, may be the twice of the eclipse time 2? eclip , which is typically ~ 3 hr. The 
timescale to release the heat in the skin depth when the particle is in Saturn's shadow (i,e., 
thermal relaxation time for eclipse) is given by 

4nrH s pC(T p - 7» _ Y(T P - 7» 

' rel " 4^v SB 6(r p 4 - r p V - «r SB (r p 4 - t; ) ^ {r ^ (2) 

where r is the radius of the particle, T p is the particle temperature, T p0 is the particle temperature 
calculated only using the Saturn's thermal flux, <x S b is the Stefan-Boltzmann constant, e is the 
thermal emissivity, and Y = yjpCK is the thermal inertia. The temperature T V) o is expected to 
be close to the temperature at equinox (when the solar elevation is zero), which occurred in 
August, 2009 (Spilker et al., 2009; Flandes et al., 2010). If t re \ «: t eclip , T p quickly drops to r p 
after particles enter Saturn's shadow. On the other hand, if t m \ » £ ec ii p , the temperature drop is 
negligible. For the typical thermal inertia values of Saturn's ring particles reported in previous 
works (~ 10 Jm~ 2 K~ 1 s" 1/2 ; Ferrari et al., 2005), t re \ ~ t ec n p , so the degree of the temperature 
drop sensitively depends on Y. 

Eq. assumes that the particle size is sufficiently larger than the skin depth (r » £ s ). If 
r «; £ s , the thermal relaxation time is given by the larger one of the thermal diffusion time 

r 2 

*diff - vl „ (3) 
K/pC 

or the time to release the internal energy of whole the hemisphere by radiation 

_ | nr 3 pC(T p - r p , ) = 1 pCr{T p - T p , ) 
tmd ~ 4nr2cr SB e(T* - T* fi ) ~ 3 ecr SB (T^ - T^) ' 

The ratio of the two time scale is roughly given as tm/tmd ~ o~sBT p r/K. Substituting typical 
values for Saturn's rings and considering a particle comparable to the skin depth (K = 10~ 3 W 
nT 2 K, T p = 90 K, r = £ s ~ 1 mm), one can find that f diff «; t md . This means that the interior of 
a small particle quickly reaches to isothermal due to heat conduction before the particle cools 
down by radiation. Thus, the temperature drop of small particles in the shadow depends on 
the product pCr instead of Y in Eq. ([2]). When r ~ £ s , t re \ depends on both pCr and Y. In 
the present paper, we first consider particles sizes sufficiently larger than the skin depth and 
estimate the thermal inertia values. Then, we make additional fits where size of small particles 
is varied. Comparison between Eqs. ([2]) and ([4]) suggests that decreasing r has an effect similar 
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to decreasing Y (i.e., particles tend to react to the input illumination instantaneously in both 
cases). We will discuss this relation in Sec. 3.3.3 after we estimate Y and r. 

Our thermal model assumes a bimodal spin distribution consisting of slow and fast rotators, 
while an actual spin distribution of a ring may be continuous. A discussion for the criterion for 
fast and slow rotators is similar to the above (Farinella et al., 1998). The time scale of the 
illumination change for a particle spin is the spin period t spin and the corresponding frequency 
is given by co = 2n/t spin . The particle is regarded as a fast rotator if ? spin «: f rel , and as a slow 
rotator if r spin » t Te \. Dynamical simulations for Saturn's rings (Ohtsuki, 2005; Morishima and 
Salo, 2006) showed that, for Saturn's rings, the spin period of the largest particles is comparable 
to the orbital period ? or b and is usually longer than the relaxation time (Y sp i n ~ t orb > t ie \), while 
the spin period of small particles is much shorter, as the spin period is roughly proportional to 
their sizes due to mutual collisions. The fraction of fast rotators / fast (the normal optical depth 
of fast rotators relative to the total normal optical depth, t), which is one of the parameters in 
our model, may be approximately given by the total cross section of particles with t sp i n < t m \ 
relative to the total cross section of all particles. Since t m \ depends on Y, / tast depends on Y, if 
a spin (and size) distribution is given. In this paper, we first estimate / fast and Y as independent 
parameters, as we did in Papers I and EL Then, we compare estimated values of / fast and Y with 
the relation predicted from the above theory using the size distribution estimated in previous 
observations (French and Nicholson, 2000) (see Sec. 4.3). 



3. Parameter fits 
3.1. Data selection 

The azimuthally scanned CIRS data for the main rings are summarized in Tables 1-4 (ob- 
servational geometries) and Figs. 1-4 (observed temperatures). Only the scans selected for [Figs. 1- 
parameter fits (see below) are shown in these tables and figures, although there are more az- 4] 
imuthal scans. All azimuthally scanned data up to Feb. 18 (DOY 49), 2007, including those [Tables 1- 
in our Tables 1, 2, and 4, are shown in Leyrat et al. (2008b). We use some additional data 4] 
obtained later: the C10 scan for the C ring (Table 1) and all the data for the Cassini division 
(Table 3). The ring spectra from the CIRS far-infrared channel, named Focal Plane 1 (FP1) 
are used. As done in all previous works using FP1 spectra (e.g., Spilker et al., 2005, 2006), 
the ring temperature is obtained by fitting a black body spectrum multiplied by a scale factor 
to the observed spectrum between 100-400 cm 1 (100-25 /mi). The scale factor includes the 
geometric filling factor, the thermal emissivity, and a scalar factor that arises when observing a 
system comprised of particles at more than one temperatures (Altobelli et al., 2007, Paper II). 
The spectral resolution for the data ranges from 0.5 to 15.5 cm 1 (Flasar et al, 2004). 

We select the data for parameter fits taking into account the following points. 



1. Data at high solar elevation are favorable (if they exist) since a large temperature drop in 
Saturn's shadow allows us to estimate the thermal inertia accurately. For the A ring and 
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the Cassini division, the eclipse time is too short if solar elevation is large and the data at 
intermediately high (~ 15°) or lower solar elevations are used. 

2. In order to resolve the degeneracy between the albedo and the fraction of fast rotators, it 
is necessary to make fits to data at low and high solar phases at the same radial location 
simultaneously. It is desirable that spatial resolutions and solar elevations of these scans 
are similar. 

3. The data for the lit face only are used for the A and B rings. The temperature of the unlit 
face of an optically thick ring is determined by how efficiently the heat is transported 
from the lit face to the unlit face and is controlled not only by the thermal inertia but also 
by the fraction of particles which are hindered from their vertical excursions through the 
midplane by mutual collisions (Paper II). Therefore, there will be a strong degeneracy 
between the thermal inertia and the fraction of bouncing particles if we include the unlit 
face data. We have a more detailed discussion in Sec. 4.2. 

4. For the A ring scans, there are azimuthal modulations not only in the absolute thermal 
brightness but also in the physical temperature due to wakes (Leyrat et al., 2008b; Ferrari 
et al., 2009). The exact reason why wakes cause temperature modulations is not well 
known, but this effect is most pronounced when the elevation angle of an observer, \B\, 
is ~ 10° as well as the brightness modulation. Therefore, we need to use data with 
\B\ <sc 10° or \B\ » 10°. For the B ring, the temperature modulation is much smaller even 
at \B\ ~ 10°, so we do not remove data by B values. 

We find that there are not many data sets available with the above restrictions. The data 
sets in four radial locations will be used (one location for each ring): 83,000 km, 105,000 km, 
120,000 km, and 129,000 km for the C ring, the B ring, the Cassini division, and the A ring, 
respectively. Since the scans CD1 and CD2 were obtained in the same day and these data have 
little overlaps in local hour angles, we merge these two scans to a single scan in parameter fits. 
The azimuthal coverages in the scans CD4, CD5, and CD6 were originally slightly larger than 
2n. We removed some of footprints in these scans in order to avoid azimuthal overlaps in a 
single scan (see Fig. 3). This is just for convenience for our fitting code, and we believe that the 
impact of removing a small fraction of data on estimated values of parameters is very small. 

The optical depth, r, of each footprint is obtained by averaging a high resolution r profile 
over the CIRS FP1 footprint (see Appendix A of Paper II) j^] The high resolution r profile from 
Voyager PPS (Esposito et al., 1983) is used for all rings and another profile from Cassini UVIS 
(Colwell et al., 2006, 2007) is also used for the A and B rings for comparison. The mean value 
of r averaged over all footprints for each scan is shown in Tables 1 to 4 as Tpps and tuvis- These 



2 In Paper II, the observational geometries used in t smoothing are those for occultations by PPS and UVIS. 
Therefore, the smoothed r is uniquely determined with a given saturnocentric radius, r p . On the other hand, in this 
work, the geometries of CIRS observations are used. Thus, the smoothed r is not always the same even if r p 's of 
the footprint centers are the same. The use of CIRS geometries improves fits particularly for the Cassini division. 
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optical depths are used in simulations. The larger value of Tuvis than r P p S is due to wakes. 
Although our model does not take into account spatial non- uniformity due to wakes, we check 
how estimated parameters vary with r for the A and B rings. 



3.2. Fitting procedures 

We use our multilayer model developed in Paper I and improved in Paper II. Our model 
solves the equation of classical radiative transfer both in thermal and visible light, taking into 
account all the heat sources for ring particles including the thermal and reflected visible fluxes 
from Saturn and the inter-particle heating and scattering. A plane-parallel approximation is 
adopted, so spatial non- uniformity due to wakes is ignored. Our model adopts a bimodal size 
distribution consisting of small fast rotators which isotropically emit thermal emissions and 
large slow rotators represented by non-spinning smooth Lambertian particles. Fast rotators are 
assumed to have a larger vertical extension than that for slow rotators. Note that the zero-volume 
filling factor of a ring assumed in classical radiative transfer, a plane -parallel approximation, and 
smooth spherical particles are probably inappropriate for Saturn's rings, and fitted parameters 
may be potentially modified in more advanced models (see also Paper II). The thermal emission 
from a ring is calculated by integrating thermal emissions from a large number of individual fast 
and slow rotators. About 3000 particles are used in a single simulation, which is roughly three 
times larger than those used in Paper II. The most important parameters are the bolometric 
Bond albedo, A v , the fraction of fast rotators, /f ast , and the thermal inertia, Y. These three 
parameters are simultaneously estimated in model fits. Since A v and /fast were well estimated 
in Paper II and ring temperatures vary smoothly and monotonically with these two parameters, 
coarse grids (the grid sizes are typically 0. 1 for A v and 0.2 for /f ast ) with a few points around 
the values estimated in Paper II are used for them. On the other hand, fine grids are used for 
T. Ring temperatures in coarse parameter grids are interpolated to those for finer grids for the 
chi- square fits shown below. 

There are some other important input parameters: the size of fast rotators r tast , the optical 
depth r, the ratio of the vertical scale height of fast rotators to that of slow rotators h r , and 
the type of the vertical motion (whether particles rebound at the midplane or not). As did in 
Paper II, we vary these parameters and check how sensitively estimated values of A v , /fast, and 
T are affected. Although there are some uncertainties in these other parameters, the base-line 
cases, for which the estimated values of parameters are marked boldly in Table 5, give the 
most reasonable results in Paper II and this workj^] The error size of simulated temperatures [Table 5] 
is estimated to be less than 0.2 K (while 0.5 K for Paper II). The physical density of 450 kg 
m" 3 and the specific heat of 760 J kg" 1 K" 1 are used as in Papers I and II. The total time for 
simulations, by which we estimated parameters shown in Table 5, was about two cpu years. 



3 In the base-line cases, we assume that h r = 3, except h t = 1 for the B ring, based on dynamical studies (Salo 
et al., 2001; Morishima and Salo, 2006), and r tast = 10 cm. This size of fast rotators is sufficiently larger than the 
skin depth, 4- The size of slow rotators is fixed to be 5 m in all cases. 
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For simulations of a single scan, we use single values of r p and r, both averaged over all 
footprints. This is valid as the dispersions in r p and r are usually small enough. Exceptions are 
scans for the Cassini division, where small changes in radii of footprint centers or spacecraft 
ranges result in large changes in smoothed r's because the r profiles, obtained in high radial 
resolution occultations, show large radial variations in this region. Therefore, we take into 
account the t variation in different footprints for scans of the Cassini division by conducting 
two simulations with low and high r's and linearly interpolating calculated spectra to the one 
for a footprint. One clear effect of the t variation is seen in the CD6 scan (Fig. 3), in which 
the temperature drops near the noon. This is caused by the mutual shadowing effect, as r near 
the noon is largest. Our model succeeds in reproducing this temperature drop (fitting curves are 
explained in Sec. 3.3). 

We make about 100 azimuthal bins and average the elevation B and longitude Cas of 
the Cassini spacecraft over the bins. These averaged values are used in simulations. Fine- 
size bins are used near the shadow lines, as azimuthal temperature variations are large around 
there. We first calculate thermal spectra in the range of 100-400 cm 4 , which is the same rage 
used in deriving observed temperatures, with the spectral resolution of 6 cm" 1 . Then, these 
spectra are smoothed over the CIRS FP1 footprints (elongated ellipses on a ring plane) with 
a gaussian weight]^] This smoothing needs to be done not only in the azimuthal direction but 
also in the radial direction. Since we use a single averaged r p value for each scan, we assume 
that ring spectra are uniform in the radial direction and vary only in the azimuthal direction in 
a footprint. This allows us to avoid conducting different simulations at slightly different radii 
and save computational time. Applying a black body fit to the smoothed spectrum, we obtain 
the ring effective (or physical) temperature. In Papers I and II, we did not make smoothing of 
spectra obtained by simulations because we analyzed thermal data which were azimuthally far 
from Saturn's shadow and azimuthal temperature variations in footprints were small there. 

Using the effective temperatures from observations and simulations, the reduced chi-square 
x\ in the three dimensional parameter space (Ay, /f ast , Y) is evaluated. We first calculate the 
reduced chi-squares for low phase X\ 0Vi and high phase xligh separately, and take the average of 
them: 

2 , A r 1 \ 1 /^obs,n — ^sim,rt(/W' ./fast* 0\ /s - . >. , C \ 

xf^iAy, /fast, H = > (for a < or crit ), (5) 

^Viow _ M ^-j* \ (T obs „ / 

7 1 V \ / ^obs /? — ^sim raO^-V, ./fast, -T) \ 

Xhi»h( A v, /fast, r) = > : — (for a > a cA ), , (6) 

■Whigh ~ M ~f\ & 0bs,H / 



4 The full-width half-maximum (FWHM) of CIRS FP1 is 2.51 mrad with a Gaussian fit (Flasar et al. 2004) and 
the diameter we assume is 5.2 mrad, which covers more than 99% of the signal. The azimuthal extension of the 
diameter projected on a ring is typically 2-3 degs but exceeds 10 degs for some data with low spatial resolutions. 
The smoothing procedure is particularly necessary for numerical data whose corresponding observational data 
have such large azimuthal extensions. 



- 11 - 



xlAw, /fast, T )=\ (*L( A V, /fast, O + ^igh( A V, /fast, O) , (7) 

where Mow and Afhigh are the total number of footprints from multiple azimuthal scans at low and 
high phases, M(= 3) is the number of parameters, n is the index for a footprint, r obSjn and r sim; „ 
are the observed and simulated temperatures for the n-th footprint, cr obs>n is the standard error 
size of the observed temperature (see below), a is the solar phase angle, and a crit is the angle 
which discriminates low and high phases. We set the error size <x obsn to be larger one of the one 
obtained from black body fits or the standard deviation of temperatures in azimuthal bins; in 
most cases, the latter is several times larger than the former (0.2-2. OK). We choose a cr i t = 83.71° 
because a slow rotator illuminated by a single heat source has stronger thermal emission than a 
fast rotator when a < a cr n and vise versa. Values of fitted parameters do not sensitively depend 
on choice of ar crit . The combination of best-fit values of (A v , /f as t, O is derived to minimize x\- 
The error bars are estimated after Press et al. (1986) as follows. The chi-square subtracted by 
its minimum value at the best-fit parameters is defined as A^. The lower and upper limits of 
a parameter is given by the maximum and minimum values of the parameter on the surface of 
A*o = 1 (68% confidence level of a parameter without regard to other parameters). When the 
upper or lower limits of parameters are out of coverages in simulations, we extrapolate x\ to 
some degree. The B ring data are split into two sets, one at high solar elevations (the scans Bl- 
B4) and another at low elevations (the scans B5- B8), and^'s are calculated separately and the 
average^ 2 is taken. The B ring parameters are also estimated for high and low solar elevation 
cases separately to see whether our model can reproduce the dependence on solar elevation. 

We first assume that all particles have the same value of T. Later, we consider a case 
with different values of Y for slow and fast rotators (r s i ow and r fast ) and parameter fits will be 
done similar to the above but in the four dimensional parameter space (M = A). The reduced 
chi-square for this case is defined as;^ 2 (A v , /fast, r s i ow , r fast ). 

In order to see the goodness of the estimated parameters, the reduced chi-square for an 
individual scan is also calculated as 

2/ a r t-\ 1 \ 1 /^obs.rc — ^sim,n(^V, /fast, O \ , r , N /ON 

X (Av, /fast, D = 77 r; > , (for each scan), (8) 

Nfoot - M j-f \ o- obs ,„ / 

where /Vf 00 t is the number of footprints in each scan shown in Tables 1-4. 



3.3. Results 

3.3.1. A case with a same thermal inertia for all particles 

We first assume that all particles have the same Y. The coverage of Y in simulations is 
2-79 Jnr 2 K-V 1/2 with twelve grid points (r,- = i{i + l)/2 + 1 Jnr 2 K-V 1/2 for i = 1 to 12). 
The estimated values of A v , /f as t, and Y for all cases are listed in Table 5 as A v ,o, /f as t,o> and 
r , respectively, and those parameters for the base-line cases are shown in Fig. 5. The best fit [Fig. 5] 
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temperature curves for the base-line cases are given by black solid curves in Figs. 1-4. The 
values of T in units of Jm _2 K _1 s _1/2 for the base-line cases are: 20.0 for the C ring, 13.0 for 
the B ring ,11.0 for the Cassini division, and 16.2 for the A ring. The values of A v ,o and /f ast; o 
are consistent with those estimated in Paper II (see Sec. 4.2 for a more detailed discussion). For 
the case of the B ring, A v ,o estimated with the high \B'\ data is lower than that with the low \B'\ 
data. This means that our model underestimates the \B'\ dependence of the ring temperature 
as compared with the observations. Nevertheless, the differences of the estimated parameter 
values for the high and low \B'\ cases remain reasonably small and the parameters estimated 
using all data simply take intermediate values between those for the high and low \B'\ cases. In 
the followings, we will mainly show the results for the high and low \B'\ cases separately, as the 
reduced chi-square values are much smaller than those for the case using all data. 

The contour of the reduced chi-square, A^, sliced on the A v vs. T plane with / fast = / fastj0 
is shown in Fig. 6. If the dispersion of observed temperatures at a given range of azimuthal [Fig. 6] 
angle is large, Axl becomes flat around the best fit parameters and it makes error sizes larger. 
The Cassini division is a such case; the large temperature dispersion is probably due to com- 
plicated radial structures whose scales are smaller than the CIRS footprint size. The dispersion 
of observed temperatures is usually small for optically thick rings due to good signal to noise 
ratio, so that the B ring parameters have small error bars. However, the actual uncertainties of 
parameters for the B ring are probably larger than at least the difference between the values for 
the high and low \B'\ cases which are larger than the error bars (particularly for A v ). 

3.3.1.1. Dependence on phase angle 

Comparisons of the observed temperatures with the modeled temperatures (black curves) 
in Figs. 1-4 show that the ring temperatures outside the shadow at both low and high phases are 
well reproduced by the model with appropriate combinations of A v ,o and /f astj0 whereas there are 
systematic deviations between the modeled and observed temperatures in the shadow depending 
on a. At low a, the shadow temperature is overestimated in the model, or T is overestimated 
(see e.g., the CI scan). On the other hand, the shadow temperature is underestimated in the 
model at high phase (e.g., the C3 scan). These facts seem to indicate that T increases with a. In 
order to confirm this trend, we make test fits using fixed A v and / fast for all scans and varying 
only T for different scans. Note that the temperature difference between inside and outside the 
shadow is controlled almost solely by F whereas a change in A v or / fast causes nearly a constant 
temperature change at all azimuthal locations. 

The estimated values of T of individual scans for the base-line cases are shown in Fig. 7 
and listed in Table 6 as r ind . The corresponding temperature curves are given by blue dashed 
curves in Figs. 1-4. As we expected, r ind increases with a for all rings. These trends are verified 
as the reduced chi-square, x 2 (shown in Figs. 1-4), for each scan is reduced with a-dependent 
T. In Fig. 7, we also plot the Y values for the C and B rings estimated by Ferrari et al. (2005). 
The range of a in their ground based observations is 4.7 - 6.2°, and their very low T values are 
consistent with our fitting lines for a vs. T. For the B ring, if we use A v and / fast estimated using 
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all data, the low \B'\ case shows lower r ind than that for the high \B'\ case at similar a (shown 
only in Table 6 not in Fig. 7), because the temperature outside the shadow is underestimated for 
the high \B\' case and overestimated for the low \B'\ case. On the other hand, the use of different 
combinations of A v and /f ast gives values of ri nd consistent each other for both cases. 

Although Ti nd depends on a as seen in Fig. 7, it is unlikely that individual particles have 
a-dependent T and their low-r faces align toward the Sun because particles are spinning and 
spin axis orientations change by collisions. A plausible interpretation for the a-dependence 
of r ind is that large slow rotators have lower T values than those for small fast rotators. The 
thermal emission from slow rotators is relatively stronger than that from fast rotators at low 
phase and vice versa. Therefore, particle groups which are dominant in the thermal emission 
may vary with phase angle and measured thermal inertias may vary as well. This possibility 
will be examined in Sec. 3.3.2. 

3.3.1.2. Dependence on optical depth and particle bouncing 

For the A and B rings, Y is almost independent of the scale height ratio, h r , but decreases 
with increasing optical depth, r, and takes larger values for the case without bouncing than 
with bouncing. The latter two dependences indicate that the thermal relaxation time for an 
optically thick ring depends not only on T of individual particles, but also on r and type of 
vertical motion, whereas the discussion for the thermal relaxation time in Sec. 2 is for an isolated 
particle. In order to help understand the latter two dependences, we examine how the ring 
physical temperature decreases in Saturn' shadow at different vertical locations, -3 < z/h < 
3, where z is the vertical position and is positive on the lit face and h is the vertical scale 
height (Fig. 8). Here, the ring physical temperature is given by the averaged temperature over [Fig. 8] 
different particles at a certain z (Eq. (27) of Paper I), and all particles are assumed to be fast 
rotators (/f ast = 1) so that the temperature is independent of the direction of emission. For 
cases with bouncing of particles at the midplane (the top and middle panels of Fig. 8), the 
energy is transported through the midplane by radiation only, so the unlit face temperature 
(z < 0) is very low and its azimuthal variation is also very small. On the lit face, the ring 
temperature increases with increasing z when illuminated by the Sun while it decreases with 
increasing z in Saturn's shadow. These trends are explained by the inter-particle shadowing and 
heating effects, respectively. The inter-particle heating cases a sort of the green house effect, 
by which an optically thicker ring can retain the heat relatively longer in the planetary shadow 
as compared with a thinner ring (compare the top and middle panels of Fig. 8). Since the 
effective thermal relaxation time of a ring increases with increasing t for a fixed T, an estimated 
T needs to decrease with increasing r in order to reproduce the observed temperature drop 
in the shadow. When particles have sinusoidal vertical motion without bouncing, the vertical 
temperature gradient is small both for the lit and unlit faces (the bottom panel of Fig. 8). In this 
case, the heat is efficiently removed from both faces, meaning a short thermal relaxation time, 
as compared with the case with bouncing. Thus, the estimated T needs to be relatively larger. 
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3.3.2. A case with different thermal inertias for slow and fast rotators 

In the last section, the thermal inertia value for all particles was assumed to be the same. 
Here, we make another fit assuming different values of the thermal inertias for slow and fast 
rotators, r slow and r fast . This means that the number of the free parameters is now four instead 
of three. The grids of r slow and r fast used in simulations are 2, 5, 9, and 15 Jm _2 K _1 s _1/ ' 2 and 10, 
20, 35, 55, 80, and 120 Jm _2 K _1 s _1/ ' 2 , respectively. The estimated values of A v fust, r s i ow , and 
r tast are shown in Fig. 5 and listed in Table 5. The corresponding temperature curves are given 
by red solid curves in Figs. 1-4. The estimated values of A v and / fast are almost the same with 
those for the case with a single thermal inertia, r . The estimated values of r slow and r fast for 
the base-line cases are in the range of 5-9 Jm~ 2 K _1 s~ 1/2 and 19-77 Jm _2 K _1 s _1/ ' 2 , respectively. 
Although the values of Ff ast are large, they are still much lower than that for solid ice with 
no porosity (~ 2600 Jm~ 2 K _1 s ~ 1/2 ). Therefore, fast rotators are likely to have surface regolith 
layers, but they may be not as fluffy as those for slow rotators. The most likely reason for 
larger r fast than r slow is that the centrifugal force for fast rotators is so strong that fluffy regolith 
particles cannot stay on fast rotators. 

Unfortunately, the error bars for r s i ow and r fast are larger than those for r , because of the 
increased freedom in adjusting temperature drops in the planetary shadow: although slow rota- 
tors are dominant emitters at low phase, the contribution from fast rotators is not unimportant 
and similarly the slow rotators' contribution is not unimportant at high phase. Figure 9 shows 
contours of Ay 2 on the r slow vs. r fast plane with the best-fit values of A v and / fast . One can [Fig. 9] 
find that t±x\ is very flat along the r fast axis and the error bars for r fast are not well determined 
for some cases. This happens particularly if / fast is small (B ring, low \B'\) or the scatter in the 
observed temperature is large (the Cassini division). However, it is obvious that r fast > r slow for 
all rings. 

Parameter fits are clearly improved by using different values of r s i ow and r tast , as the lowest 
x\ (with the best-fit parameters) is lower than the lowest x\ f° r ai l rings. Indeed, the chi- 
square for an individual scan, x 2 , is also reduced for almost all scans for all rings. Although 
the a-dependent systematic deviations between the modeled and observed temperatures in the 
shadow are much reduced, they are not completely removed, as x 2 is still larger than that for the 
a-dependent Y case. This is probably because our model assumes a bimodal size distribution 
whereas an actual ring has a continuous size distribution. The distribution of the thermal inertia 
is likely to be continuous as well. Therefore, Y for the largest particles may be even smaller than 
r s iow> which is the thermal inertia averaged over a certain size range of large particles. Similarly, 
T for the smallest particles may be larger than r fast . This effect needs to be investigated in future 
works. 



3.3.3. Lower limit of the size of fast rotators 

As discussed in Sec. 2, fitting results shown in Sees. 3.3.1 and 3.3.2 are not altered even if 
we use different particle sizes as long as the particle size is sufficiently larger than the thermal 
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skin depth. If we use particle sizes as small as or even smaller than the skin depth, reducing 
particle size causes an effect similar to reducing thermal inertia. Eqs. Q and (|4]) suggest the 
following relation: 

3r 

r«— -=. (9) 
pC yto 

Since the use of a too small r makes fits worse, the lower limit of the size of particles can be 
constrained. Here we make additional parameter fits varying size of fast rotators, r fast . The num- 
ber of parameters to be fitted is now five instead of four in Sec. 3.3.2. We define the lower limit 
size of fast rotators, rf ast m i n , as the lowest rf ast giving A^ 2 = 1 in the five dimensional parameter 
spacej^] This is similar to the definitions of the lower limits of other parameters estimated in 
previous sections. Since a comprehensive coverage of the five dimensional parameter space is 
computationally so intense that the following effective method to estimate r fast min is adopted. 
We first fix all four parameters (Ay, /fast, r s i ow , and r fast ) estimated in Sec. 3.3.2 and search for a 
size of fast rotators, r 2 , which gives A^ 2 ~ 1. Then, we make parameter fits varying four other 
parameters at r tast = r 2 and 0.5r 2 . When other parameters vary, is usually smaller than 1 for 
ffast = r 2- This means that r fast min should be at least smaller than r%. Linearly interpolating Ax] 
for '"fast = r 2 and 0.5r 2 , we estimate r fast min , where Ax\ = 1. 

The values of ^ for rf ast ~ r 2 and other estimated parameters are listed in Table 5. We find 
that r fastimin ranges from 1-10 mm. The temperature curves for r fast ~ rf astimin are shown as purple 
dotted curves in Figs. 1-4. Figure 10 compares r fastmin with the lower limit of r fast , r fast;mill , 
for the base-line cases (r fast = 10 cm) (note if r fast = 33.5+jgg Jm" 2 K _1 s" 1/2 , r fastjmin is 18.0 
Jm _2 K _1 s _1/2 ). For the B ring at low \B'\ and the Cassini division, we were not able to estimate [Fig. 10] 
Ffast,min- For other three rings, for which the values of Ff astim i n are well estimated, the relations 
between r fast , min and r fast , min are well represented by Eq. ^ with r = r fastimin and Y = r fast;min . 
There, the cycle of the illumination change, 2n/a>, needs to be specified to determine the thermal 
skin depth. We find that a cycle which gives a reasonable fit is somewhere between the twice 
of the eclipse time 2f eclip (see Tables 1-4 for t ec n p ) and the orbital period f orb J^] We think that this 
is reasonable, because temperatures of ring particles recover almost immediately after they exit 
Saturn's shadow, but not completely, and their temperatures gradually increase until they enter 
Saturn's shadow again. 



Correction factors (Mow -4) /(Mow -5) and (Mngh-4)/(Migh-5) are ignored, as Mow andMigh are sufficiently 
large. 

6 The orbital periods are 6.78 hr for the C ring (83,000 km), 9.65 hr for the B ring (105,000 km), 1 1.79 hr for 
the Cassini division (120,000 km), and 13.14 hr for the A ring (129,000 km). 
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4. Discussion 
4.1. Interpretations 

The estimation of r fastim i n in Sec. 3.3.3 means that particles with sizes equal to or less than 
7fast,min should not occupy more than / fast of the total cross section. This condition seems to be 
satisfied for all rings, as the lower cutoff r min for the continuous size distribution estimated in 
French and Nicholson (2000) (Table 7) is comparable to or larger than r fast>min . Therefore, the [Table. 7] 
quantity constrained by the eclipse cooling is the thermal inertia, not the size of ring particles. 

In Sec. 3.3.2, we found r fast > r s i ow for all the rings. The most likely reason for this 
is that small fast rotators cannot hold very fluffy regolith layers because of their fast spins. 
Dynamical studies (Salo, 1987; Ohtsuki, 2005; Morishima and Salo, 2006) show that the spin 
period is about the orbital period for the largest particles and is roughly proportional to the 
particle size for an extended size distribution. The centrifugal force is comparable to the self- 
gravity on a particle with a synchronous rotation near the Roche limit, around which Saturn's 
main rings exist. Therefore, the centrifugal force well exceeds the self-gravity on surfaces of 
small particles. Since r fast is likely to be still much smaller than the thermal inertia of the solid 
ice with no porosity, some regolith particles probably stick to small fast rotators due to the 
surface energy (Choksi et al., 1993, Albers et al., 2006). The regolith layers of fast rotators may 
also be compacted by mutual impacts whose velocities are larger than the velocity dispersion 
of large particles (Salo, 1992; Morishima and Salo, 2006). Alternatively, small particles might 
be regolith free. In this case, as suggested by Kouchi et al. (1992), microcracks would be 
responsible for r fast . In any case, it is expected that the thermal inertia increases with decreasing 
size of a ring particle. A similar size dependence of the thermal inertia is also found in asteroids 
(Delbo and Tanga, 2009). However, the mechanism causing the size dependence is probably 
different from that for Saturn's ring particles; their interpretation is that large asteroids (> 100 
km) are primordial and collisional debris build up thick regolith layers on them whereas small 
bodies are fragments in catastrophic impacts between large bodies. 

In the interpretation of the variations of r s i ow and r fast in different rings, we need to consider 
(1) the distribution of T(r) as a function of particle size r and (2) the ratio of the smallest particle 
size r m i n to the largest particle size r max for actual rings with continuous size distribution. As 
we discussed in Sec. 3.3.2, the thermal inertia of the largest particles, T(r = r max ), may be close 
to r slow or slightly less. Therefore, T(r = r max ) seems to be similar for all the rings, except the 
Cassini division has a lower T(r = r max ). Because F is likely to increase with decreasing r, as 
discussed above, r fast is expected to increase with decreasing r min /r max . On the other hand, r fast 
for the A ring is larger than that for the C and B rings, even though r m i n /r max (A ring) = r m i n /r max 
(B ring) > r min /r max (C ring) (Table 7). This means that with a same r/r max (< 1), an A ring 
particle has the largest T. Thus, the size dependence of T for the A ring is steeper than those for 
the B and C rings, whereas r fast for the C ring larger than that for the B ring may be explained 
by the smaller r min /r max for the C ring. The large r fast for the A ring is probably due to wakes 
(Salo, 1995; Colwell et al., 2006; Hedman et al., 2007). Dynamical simulations (Morishima and 
Salo, 2006) including the self-gravity show that the spin frequency increases with increasing r p 
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for a fixed dynamical optical depth. The self-gravity also steepens size dependence of the spin 
frequency. The collision velocity between large particles are usually slow as they are in the same 
wake, while small particles floating in inter-wake spaces suffer collisions with impact velocities 
as fast as the velocity dispersion of wakes, which is much larger than the escape velocity of 
individual particles. Therefore, for small particles with a certain r/r max («: 1), both the spin 
frequency and the collision velocity are largest in the A ring. These facts may explain large r fast 
for the A ring, while r slow for the A ring is similar to those for the C and B rings. It should be 
noted again that the spatial non- uniformity due to wakes is not directly taken into account in 
our model. The discussion here was based on the estimated values of r fast for the standard cases 
and readers are cautioned that the uncertainty in r fast is very large. 

The Cassini division has the smallest r slow among the Saturn's main rings. This smallest 
value seems to be real, as r ind in Fig. 7 at a given a always shows the smallest value, except at 
very large a (i.e., the CD3 scan). The very small value of r s i ow may indicate that there are some 
dust sources around this region and ring particles are fluffily coated by these grains. There are 
at least two dusty ringlets in the Cassini division (Horanyi et al., 2009): one in the Huygens 
gap (117,490 km) and another in the Laplace gap (119,940 km) called the Charming ringlet 
(Hedman et al., 2010). In particular, the location of the latter one is very close to the centers 
of the CIRS footprints of azimuthal scans. These dusty ringlets themselves are very tenuous 
and unimportant for the direct thermal emission. However, their existence indicates injection 
of small grains into this region by some unknown mechanisms, as they are likely to be very 
young (less than a few thousand years old; Hedman et al., 2010). The small lower-cutoff size 
for the Cassini division found in French and Nicholson (2000) (Table 7) also seems consistent 
with a relatively large population of dust. They also inferred that there is a substantial variation 
in the particle size distribution across this region, with relatively more small particles in the 
inner than in the outer Cassini division (reported also in Colwell et al. (2009)). This probably 
indicates that the dust sources are limited to the inner Cassini division. The value of r ind for 
the CD3 scan which looks remarkably too large as compared with those for other scans seems 
to be explained by the fact that the saturnocentric distance of the CD3 scan is slightly larger 
than those for other scans (by ~ 2000 km at shadow lines). Howett et al. (2010) found that 
the thermal inertias of saturnian satellites from Rhea inward have relatively lower values (~ 10 
Jm~ 2 K _1 s~ 1/2 ) than those for Iapetus and Phoebe (~ 20 Jm _2 K _1 s ~ 1/2 ), which are located farther 
away from Saturn. They indicated that this is probably related to coating by E-ring material 
ejected from Enceladus' plumes. Similar things may happen also in Saturn's rings although 
dust sources are different. 

The large value of x\ f° r the B ring in the case using all data comes from the fact that 
our model underestimates \B'\ dependence of the ring temperature, as we stated in Sec. 3.3.1. 
By separately estimating parameters for the B ring at the high and low \B'\ cases, the values of 
x\ are remarkably reduced. However, x\ for the high \B'\ case is still much larger than unity. 
The phase angles for the scans B2, B3, and B4 are low. In Fig. 2, with decreasing phase angle 
from ~ 40° (the B4 scan) to ~ 10° (the B2 scan), the observed temperature outside the shadow 
increases by ~ 4K, whereas the modeled temperature increases by only ~ IK. Therefore, our 
model underestimates not only \B'\ dependence but also a dependence at low phase. The steep a 
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dependence observed is probably due to the thermal opposition surge caused by packing of ring 
particles into thin layers (Altobelli et al., 2009). The effect of particle packing is also known to 
enhance \B' \ dependence of photometric brightness for dense rings (Salo and Karjalainen, 2003), 
but it is uncertain if it is also the case for the thermal emission. For more accurate estimations 
of particle properties for optically thick rings such as the B and A rings, this packing effect and 
the effect due to wakes need to be taken into account in future models. Nevertheless, we believe 
that parameters estimated in the present work are not so different from actual values, because 
the deviations between observed and modeled temperatures remain reasonably small. 



4.2. Toward estimations of radial variation of thermal inertia 

We have estimated the thermal inertia values for only four radial points in this study (one 
point for each ring), and it is desirable to examine the radial variation of the thermal inertia 
in each ring. Unfortunately, azimuthal scans were taken at only a limited number of radial 
points by Cassini CIRS and it is probably difficult to cover a sufficiently large number of radial 
points even in the extended mission. However, there are some quasi-radial scans along the 
shadow lines and in Saturn's shadow. Even if we use reasonable sets of shadow scans, precise 
estimation of T is not simple in the optically thick rings (the A and B rings) because an estimated 
T depends on types of vertical motion of ring particles and adopted r (Table 5). 

The degeneracy between Y and the type of vertical motion can be resolved if we combine 
T estimated from the temperature difference between the lit and unlit faces (Paper II) and F 
estimated from the eclipse cooling (this work). Figure 1 1 shows comparisons of parameters 
estimated in Paper II (solid lines) and this work (diamonds) for the case with t pps . Here we [Fig. 11] 
use r from this study because Paper II did not consider different values of T for slow and fast 
rotators. For the B and A rings, estimated parameters for the bouncing model are shown with 
black lines. Very good agreements between both studies are found if we assume the bouncing 
model for the B ring and the non-bouncing model for the other rings. In fact, with the non- 
bouncing model, T and A v for the A ring from this study are slightly higher and lower than those 
from Paper II. These deviations are fixed if we assume a small fraction of particles rebound at 
the midplane, because both Y from Paper II and A v from this study increase with increasing 
fraction of bouncing particles. 

There is still a problem in fitted parameters due to their dependences on t. The uncertainty 
in r for the A and B rings comes from wakes (Colwell et al., 2006, 2007), so ultimately more 
sophisticated models, which take into account wakes, are necessary. However, the uncertainty 
in estimated values of Y due to the uncertainty of r seems very small, as we will see below. 
Figure 12 shows Y for the mid B ring (105,000 km) estimated from this study (the mean values 
of r for high B' and low B' cases) and from Paper II (r Pn ) as a function of the fraction of 
bouncing particles ft, for two different cases of optical depths, Tpps and Tuvis- Here, f b = and [Fig. 12] 
1 represent the non-bouncing and bouncing models, respectively, and a linear dependence of 
T on f b is assumed between f b = and 1 . In this work, we did not estimate r for the case 
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with fb = and t UV is, and we assume that the ratio r (Tpp S )/ro(T UV is) is independent of f b 
as well as the ratio of the error sizes. The values of Y and f b on the crossing points of the 
solid and dashed lines are the predicted values. For the case with Tp PS , the crossing point has 
T ^ 10-13 JnT 2 Kr 1 s~ 1,/2 and / b - 1- On the other hand, for the case with tuvis, the crossing 
point has T = 14 Jm^K^'s" 1 ^ and fb = 0.21. Since we used the possible minimum value, 40 
Jnr 2 K _1 s~ 1/2 , for r PII (Tuvis) at/t, = 1, T and f b at the cross point will be respectively higher and 
lower than those shown in Fig. 12. Nevertheless, Y at the crossing point should not be larger 
than r (T UV is) at fb = (~15 Jm ^K^s ^ 2 ). Therefore, even though we cannot resolve the 
degeneracy between r and fb, similar values of Y seem to be obtained regardless of the choice 
of r. This discussion needs to be extended to the case with different values of r s i ow and r fast in 
future works. 



4.3. Criteria of slow and fast rotators 

We have treated the thermal inertia Y and the fraction of fast rotators / fast as independent 
parameters in data fitting. However, they are not independent of each other if the spin (and size) 
distribution is given. Here, we discuss the relation between Y and /f ast . A similar discussion is 
found in Sec. 5.2 of Paper I. 

We adopt the size distribution estimated in French and Nicholson (2000) (Table 7). The 
number density per unit size n(r) is given by a power-law as n(r) oc r"' y (r min < r < r max ), where 
r mm and r max are the minimum and maximum size of particles, respectively. They estimated the 
size distribution using diffracted stellar light, and their results are consistent with those from 
radio occultation (Zebker et al., 1985; Cuzzi et al., 2009). From dynamical studies (Ohtsuki, 
2005, 2006; Morishima and Salo, 2006), the particle spin period is approximately given as 

r 

h V m(r) = Ci t orh , (10) 

''max 

where C\ is the numerical constant and depends on r p and r; we take Ci = 2.0, 1.0, 2.0, and 0.5 
for the C ring, the B ring, the Cassini division, and the A ring from Fig. 1 1 of Morishima and 
Salo (2006). The critical particle size, which discriminates slow and fast rotators, is defined as 

^spin^crit) — ^rel> (11) 

where f re i is the thermal relaxation time given by Eq. ([2]) with a> = 2n/t sp i n . The critical size r cr i t 
is given as 
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(12) 



p.i 

Here, Y is the thermal inertia for particles with r = r crit , which should be somewhere between 
r slow and r fast . In a comparison below, we assume that T = r as a nominal value. The tem- 
perature, r p o is supposed to be derived from the flux to the night-side hemisphere of a particle, 
and we take its value from the equinox data assuming the night-side hemisphere is illuminated 
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only by Saturn: 65 K for the C ring and 48 K for other rings (Spilker et al., 2009). If we simply 
regard a particle with r < r crit as a fast rotator and a particle with r > r crit as a slow rotator, the 
fraction of fast rotators in cross section / fast is given as 

r _ (^crit/^max) ^ ('"min/^'max) ^ fl'lx 

/fast " 1 - (r ■ Ir ) 3 ~<? " U } 

1 V mm / ' max/ 



Eqs.(12) and (13) indicate that / fast increases with increasing Y and deceasing T p and t orb , for 



given r^/r^ and <?. 



Figure 13 shows the relation between Y and / fast derived from Eqs. ( |T0| ) - ( [13] ) with r min /r max 
and g from French and Nicholson (2000). The estimated values, r and /f as t,o> from our parame- [Fig. 13] 
ter fits (Table 5) are also plotted. We find good agreements between values from our parameter 
fits and the theoretical prediction for all rings. In Paper II, we discussed that our model prob- 
ably underestimates ff ast due to the lack of the effects of packing of ring particles and surface 
roughness. The good agreements seen in Fig. 13 may indicate that these effects do not largely 
change estimated values of parameters such as /f ast and Y. For the A and B rings, we may need 



to use a larger T pfl due to the mutual heating. In this case, / fast predicted from Eq. (13) will 



be lower than those shown in Fig. 13 and the agreements will be even better unless T p0 is very 



close to T p . 



5. Conclusions 

In the present paper, we have estimated the thermal inertia values of Saturn's main rings 
by applying our thermal model to azimuthally scanned spectra taken by the Cassini CIRS in- 
strument. Model fits show the thermal inertia of ring particles to be: 16^°, 13^, 20^°, and 
ll*g 8 Jm" 2 K" 1 s" 1/2 for the A, B, and C rings and the Cassini division, respectively. However, 
there are systematic deviations between modeled and observed temperatures in Saturn's shadow 
depending on solar phase angle. Test fits, which allow individual scans taken at different so- 
lar phase angles to have different thermal inertia values, show that the apparent thermal inertia 
increases with solar phase angle. This dependence is likely to be explained, if large slowly 
spinning particles have lower thermal inertia values than those for small fast spinning particles, 
because the thermal emission of the slow rotators is relatively stronger than that of fast rotators 
at low phase and vise versa. We made additional parameter fits, which assume that slow and 
fast rotators have different thermal inertia values. The estimated thermal inertia values of slow 
rotators are 8^ 9 , 8^°, 9^ 3 , and 5^ 7 Jm" 2 K" 1 s" 1/2 for the A, B, and C rings and the Cassini 
division, respectively, and those of fast rotators are 27^ 7 3 , 34^, 55^j 5 Jm _2 K _1 s" 1/2 . 

The values for fast rotators are still much smaller than those for the solid ice with no porosity. 
Thus, fast rotators are likely to have surface regolith layers, but these may be not as fluffy as 
those for slow rotators, probably because the capability of holding regolith particles is limited 
for fast rotators due to the centrifugal force. The large thermal inertia of fast rotators for the 
A ring is probably due to wakes, which largely enhance spin and collision velocities of small 
particles. The low thermal inertia of the Cassini division probably indicates that there are some 
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dust sources around this region and ring particles are fluffily coated by dust. Other additional 
parameter fits, in which size of fast rotators is varied, indicate that particles less than ~ 1 cm 
should not occupy more than a half of the cross section for the A, B, and C rings. This is 
consistent with previous estimations of the ring particle size distribution. 



Acknowledgments 

We thank Cedric Leyrat and an anonymous reviewer for their comments, which improved 
our manuscript. We are grateful for the support by the Cassini project and the NASA's OPR 
and PGG Programs. K.O. is also grateful for the support by JSPS KAKENHI (22340125). 
We thank Nicholas Altobelli and Stu Pilorz for developing the CIRS database, Scott Edington, 
Shawn Brooks and Mark Showalter for designing the CIRS ring observations. R.M. thanks the 
UVIS ring team in LASP for fruitful discussions. Numerical simulations were carried out with 
the supercomputers, Nebula and Galaxy, at JPL. 



REFERENCES 

Albers, N., Spahn, E, 2006. The influence of particle adhesion on the stability of agglomerates 
in Saturn's rings. Icarus 181, 292-301. 

Altobelli, N., Spilker, L.J., Pilorz, S., Brooks, S., Edgington, S., Wallis, B., Flasar, M, 2007. 
C ring fine structures revealed in the thermal infrared. Icarus 191, 691-701. 

Altobelli, N., Spilker, L.J., Leyrat, C, Pilorz, S., 2008. Thermal observations of Saturn's main 
rings by Cassini CIRS: Phase, emission and solar elevation dependence. Planet. Space. 
Sci. 56, 134-146. 

Altobelli, N., Spilker, L.J., Pilorz, S., Leyrat, C, Edgington, S., Wallis, B., Flandes, A., 2009. 
Thermal phase curves observed in Saturn's main rings by Cassini-CIRS: Detection of an 
opposition effect? Geophys. Res. Lett. 36, L10105. 

Aumann, H.H., Kieffer, H.H., 1973. Determination of particle sizes in Saturn's rings from 
their eclipse cooling and heating curves. Astrophys. J. 186, 305-311. 

Bradley, E.T Colwell, J.E., Esposito, L.W, Cuzzi, J.N., Tollerud, H., Chambers, L., 2010. Far 
ultraviolet spectral properties of Saturn's rings from Cassini UVIS. Icarus 206, 458^4-66. 

Choksi, A., Tielens, A.G.G.M., Hollenbach, D., 1993 Dust coagulation. Astrophys. J. 407, 
806-819. 

Colwell, J.E., Esposito, L.W., Sremcevic, M., 2006. Self-gravity wakes in Saturn's A ring 
measured by stellar occultations from Cassini. Geophys. Res. Lett. 33, L07201. 



-22- 



Colwell, J.E., Esposito, L.W., Sremcevic, M., Stewart, G.R., McClintock, W.E., 2007. Self- 
gravity wakes and radial structure of Saturn's B ring. Icarus, 190, 127-144. 

Colwell, J.E., Cooney, J.H., Esposito, L.W., Sremcevic, M., 2009. Density waves in Cassini 
UVIS stellar occultations: 1. The Cassini Division. Icarus 200, 574-580. 

Cuzzi, J.N., Clark, K., Filacchione, G., French, R., Jhonson, R., Marouf, E., Spilker, L., 2009. 
Ring particle composition and size distribution. In: Dougherty, M.K., Esposito, L.W., 
Krimigis, S.M. (Eds.), Saturn from Cassini-Huygens. Springer, Berlin, pp. 459-509. 

D'Aversa, E., Bellucci, G., Nicholson, P.D., Hedman, M.M., Brown, R.H., Showalter, M.R., 
Altieri, F, Carrozzo, EG., Filacchione, G., Tosi, F, 2010. The spectrum of a Saturn ring 
spoke from Cassini/VIMS. Geophys. Res. L. 37, L01203. 

Delbo, M., Tanga, R, 2009. Thermal inertia of main belt asteroids smaller than 100 km from 
IRAS data. Planetary and Space Sci. 57, 259-265. 

Dones, L., Cuzzi, J.N., Showalter, M.R., 1993. Voyager photometry of Saturn's A ring. Icarus 
105, 184-215. 

Esposito, L.W., O'Callaghan, M., Simmons, K.E., Hord, C.W., West, R.A., Lane, A.L., Pom- 
phrey, R.B., Coffeen, D.L., Sato, M., 1983. Voyager photopolarimeter steller occultation 
of Saturn's rings. J. Geophys. Res. 88, 8643-8649. 

Farinella, P., Vokorouhlicky, D., Hartmann, W.K., 1998. Meteorite delivery via Yarkovsky 
orbital drift. Icarus 132, 378-387. 

Ferrari, C, Leyrat, C, 2006. Thermal emission of spherical spinning ring particles: The 
standard model. Astron. Astrophys. 447, 745-760. 

Ferrari, C, Galdemard, P., Lagage, P.O., Pantin, E., Quoirin, C, 2005. Imaging Saturn's rings 
with CAMIRAS: thermal inertia of B and C rings. Astron. Astrophys. 441, 379-389. 

Ferrari, C, Brooks, S., Edgington, S., Leyrat, C, Pilorz, S., Spilker, L., 2009. Structure of 
self-gravity wakes in Saturn's A ring as measured by Cassini CIRS. Icarus 199, 145-153. 

Flandes, A., Spilker, L.J., Morishima, R., Pilorz, S., Leyrat, C, Altobelli, N., Brooks, S.M., 
Edgington, S.G., 2010. Brightness of Saturn's rings with decreasing solar elevation. 
Planet and Space Science 58, 1758-1765. 

Flasar, F.M., and 44 colleagues, 2004. Exploring the Saturn system in the thermal infrared: 
The composite infrared spectrometer. Space Science Rev. 115, 169-297. 

Flasar, F.M., and 45 colleagues, 2005. Temperatures, winds, and composition in the Saturnian 
system. Science 307, 1247-1251. 

French, R.G., Nicholson, P.D., 2000. Saturn's rings II. Particle sizes inferred from stellar 
occultation data. Icarus 145, 502-523. 



-23- 



Froidevaux, L., 1981. Saturn's rings: Infrared brightness variation with solar elevation. Icarus 
46, 4-17. 

Froidevaux, L., Matthews, K., Neugebauer, G., 1981. Thermal response of Saturn's ring par- 
ticles during and after eclipse. Icarus 46, 18-26. 

Hedman, M.M., Nicholson, P.D., Salo, H., Wallis, B.D. Buratti, B.J., Baines, K.H., Brown, 
R.H., Clark, R.N., 2007. Self-gravity wake structures in Saturn's A ring revealed by 
Cassini VIMS. Astron. J. 133, 2624-2629. 

Hedman, M.M., Burt, J. A., Burns, J.A., Tiscareno, M.S., 2010. The shape and dynamics of a 
heliotropic dusty ringlet in the Cassini division. Icarus 210, 284-297. 

Horanyi, M., Burns, J.A., Hedman, M.M., Jones, G.H., Kempf, S., 2009. Diffuse rings. In: 
Dougherty, M.K., Esposito, L.W., Krimigis, S.M. (Eds.), Saturn from Cassini-Huygens. 
Springer, Berlin, pp. 511-536. 

Howett, C.J.A., Spencer, J.R., Peral, J., Segura, M., 2010. Thermal inertia and bolometric 
Bond albedo values for Mimas, Enceladus, Tethys, Dione, Rhea, and Iapetus as derived 
from Cassini/CIRS measurements. Icarus 206, 573-593. 

Kouchi, A., Greenberg, J.M., Yamamoto, T, Mukai, T, 1992. Extremely low thermal conduc- 
tivity of amorphous ice: relevance to comet evolution. Astrophys. J. 388, L73-L76. 

Leyrat, C, Ferrari, C, Charnoz, S., Decriem, J., Spilker, L., Pilorz, S., 2008a. Spinning 
particles in Saturn's C ring from mid-infrared observations: Pre-Cassini results. Icarus 
196,625-641. 

Leyrat, C, Spilker, L, J., Altobelli, N., Pilorz, S., Ferrari, C, 2008b. Infrared observations of 
Saturn's rings by Cassini CIRS : Phase angle and local time dependence Planet. Space. 
Sci. 56, 117-133. 

Marouf, E.A., Tyler, G.L., Zebker, H.A., Simpson, R.A., Eshleman, V.R., 1983. Particle size 
distributions in Saturn's rings from Voyager 1 radio occultation. Icarus 54, 189-211. 

Morishima, R., Salo, H., 2006. Simulations of dense planetary rings IV: Rotating self-gravitating 
particles with size distribution. Icarus 181, 272-291. 

Morishima, R., Salo, H., Ohtsuki, K., 2009 (Paper I). A multilayer model for thermal infrared 
emission of Saturn's rings: Basic formulation and implications for Earth-based observa- 
tions. Icarus 201, 634-654. 

Morishima, R., Spilker, L., Salo, H, Ohtsuki, K., Altobelli, N., Pilorz, S., 2010 (Paper II). A 
multilayer model for thermal infrared emission of Saturn's rings. II: Albedo, spins, and 
vertical mixing of ring particles inferred from Cassini CIRS. Icarus 210, 330-345. 

Morrison, D., 1974. Infrared radiometry of the rings of Saturn. Icarus 22, 57-65. 



-24- 



Nicholson, P.D., and 15 coauthors, 2008. A close look at Saturn's rings with Cassini VIMS. 
Icarus 193, 182-212. 

Ohtsuki, K., 2005. Rotation rates of particles in Saturn's rings. Astrophys. J. 626, L61-L64. 

Ohtsuki, K., 2006. Rotation rate and velocity dispersion of planetary ring particles with size 
distribution II. Numerical simulation for gravitating particles. Icarus 183, 384-395. 

Poulet, E, Cruikshank, D.P., Cuzzi, J.N., Roush, T.L., French, R.G., 2003. Composition of 
Saturn's rings A, B, and C from high resolution near-infrared spectroscopic observations. 
Astron. Astrophys. 412, 305-316. 

Press, W.H., Teukolsky, S.A., Vetterling, W.T., Flannery, B.P, 1986. Numerical Recipes. 
Cambridge Univ. Press, Cambridge, UK. 

Salo, H., 1987. Numerical simulations of collisions between rotating particles. Icarus 70, 
37-51. 

Salo, H., 1992. Numerical simulations of dense collisional systems. II - Extended distribution 
of particle sizes. Icarus 96, 85-106. 

Salo, H., 1995. Simulations of dense planetary rings. III. Self-gravitating identical particles. 
Icarus 117,287-312. 

Salo, H., Karjalainen, R., 2003. Photometric modeling of Saturn's rings I. Monte Carlo method 
and the effect of nonzero volume filling factor. Icarus 164, 428-460. 

Salo, H., Schmidt, J., Spahn, E, 2001. Viscous overstability in Saturn's B ring I. Direct simu- 
lations and measurement of transport coefficients. Icarus 153, 295-315. 

Shoshany, Y., Prialnik, D., Podolak, M., 2002. Monte Carlo modeling of the thermal conduc- 
tivity of porous cometary ice. Icarus 157, 219-227. 

Spilker, L.J., Pilorz, S.H., Edgington, S.G., Wallis, B.D., Brooks, S.M., Pearl, J.C., Flasar, 
EM., 2005. Cassini CIRS observations of a roll-off in Saturn ring spectra at submillimeter 
wavelengths. Earth, Moon, and Planets 96, 149-163. 

Spilker, L.J., and 11 colleagues, 2006. Cassini thermal observations of Saturn's main rings: 
Implications for particle rotation and vertical mixing. Planet. Space Sci. 54, 1 167-1 176. 

Spilker, L.J., Flandes, A., Morishima, R., Altobelli, N., Leyrat, C, Pilorz, S., Ferrari, C.C., 
Edgington, S.G., Brooks, S.M., 2009. Saturn ring temperatures at equinox with Cassini 
CIRS. American Geophysical Union, Fall Meeting 2009, abstract #P51B-1 131. 

Stuart, S.J., Stewart, G.R., Lewis, M.C., Colwell, J.E., Sremcevic, M., 2010. Estimating the 
masses of Saturn's A and B rings from high-optical depth TV-body simulations and steller 
occultations. Icarus 206, 431-445. 



-25- 



Zebker, H.A., Marouf, E.A., Tyler, G.L., 1985. Saturn's rings - Particle size distributions for 
thin layer model. Icarus 64, 531-548. 



-26- 



Scan No. 


Date 


Afoot 


«o 


B' 


B 


a 


<P P 


r v 


<f>Cas 


tas 


TPPS 


^eclip 








(AU) 


(deg) 


(deg) 


(deg) 


(deg) 


(1,000 km) 


(deg) 


(''Sat) 




(hr) 


CI 


05-156 


1610 


9.076 


-21.48 


-18.73 


41.38 


261.10 (co) 


83.10 


143.00 


26.30 


0.108 


1.56 












-20.02 


34.67 


17.25 


82.81 


135.78 


24.46 






C2 


05-176 


1979 


9.079 


-21.28 


-20.20 


30.51 


276.91 (co) 


83.11 


157.85 


19.19 


0.111 


1.56 












-22.30 


20.64 


27.82 


82.80 


147.23 


17.08 






C3 


05-178 


1711 


9.079 


-21.25 


9.59 


125.27 


286.54 (cl) 


83.59 


59.99 


11.18 


0.107 


1.57 












9.55 


121.33 


1.69 


82.39 


55.75 


10.56 






C4 


05-194 


1764 


9.081 


-21.09 


-19.88 


28.25 


22.05 (cl) 


83.08 


151.57 


21.31 


0.110 


1.57 












-21.76 


26.49 


100.21 


82.81 


149.83 


19.03 






C5 


06-300 


1126 


9.164 


-15.29 


-22.66 


112.84 


283.54 (cl) 


83.76 


60.10 


12.92 


0.111 


1.66 












-28.19 


107.66 


1.37 


82.90 


56.61 


10.88 






C6 


06-337 


2482 


9.172 


-14.77 


60.00 


113.76 


270.33 (cl) 


83.58 


289.52 


12.77 


0.111 


1.67 












50.03 


105.73 


89.70 


82.91 


276.30 


9.98 






C7 


07-016 


1041 


9.181 


-14.16 


-48.86 


45.17 


297.13 (d) 


83.11 


171.86 


13.40 


0.112 


1.67 












-53.71 


37.09 


119.60 


82.42 


146.87 


11.96 






C8 


07-048 


1132 


9.188 


-13.70 


-25.66 


117.50 


281.25 (cl) 


84.54 


51.69 


18.54 


0.114 


1.68 












-29.92 


114.26 


2.20 


82.66 


49.87 


16.82 






C9 


07-049 


2188 


9.189 


-13.69 


-52.95 


80.11 


92.11 (co) 


83.92 


100.70 


16.47 


0.114 


1.68 












-59.23 


72.70 


301.52 


82.93 


86.40 


15.09 






CIO 


07-082 


982 


9.196 


-13.22 


-50.18 


75.22 


346.88 (cl) 


84.01 


101.22 


9.97 


0.113 


1.69 












-60.14 


72.80 


74.78 


81.49 


97.32 


9.10 







Table 1: Geometry data of azimuthal scans for the C ring. Date is year and day of year, Nf oot 
is the number of footprints (or spectra), a Q is the heliocentric distance, B' is the solar elevation 
angle, B is the elevation angle of the Cassini spacecraft, a is the solar phase angle, <7) p is the local 
hour angle (or the longitude) of ring particles around Saturn with the origin at the midnight, r p 
is the saturnocentric distance, Ca .s is the longitude of the spacecraft around the footprint, r Cas 



is the range of the spacecraft (r Sat is the radius of Saturn), t pps is the mean optical depth from 
observations taken by Voyager PPS (see the text), and t ec u p is the mean time of the eclipse. For 
B, a, r p , <7) Cas , and r Cas , the upper and lower limits of each parameter are shown in the upper 
and lower rows. For (p v , the values of initial and last foot prints during the scans are shown in 
the upper and lower rows. The characters in the parenthesis in the upper row of <p v represents 
the azimuthal direction of the scan: "co" is for a counter-clockwise scan whereas "cl" for a 
clockwise scan. 
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Table 2: Same as Table 1 but for the B ring. Here the mean optical depth from observations 



taken by Cassini UVIS, r UV is, is also shown, in addition to r PP s. 
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Table 3: Same as Table 1 but for the Cassini division. 
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Table 4: Same as Table 2 but for the A ring. 
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9.63 


6 - 8 -oi ( 6 - 3 -oD 


C4 


26.95 


11 9+ 36 

"''-2.5 


B6 


31.59 


10.0+" (7.0+JJ) 


CI 


37.68 


15 7 +5 - 9 


B3 


32.62 


10.8+j * (H- 5 -!') 


C7 


39.23 




B4 


40.05 


12.2+f i (11.9+f i) 


CIO 


73.56 


17.8^ 


B8 


83.37 


12.9^(7.1^) 


C9 


74.62 


16 6 +12 - 9 


Bl 


104.17 


16.5 +2 ;« (21.1+3 7) 


C6 


106.23 


94 <;+5(5-3 
-12.2 


B5 


104.71 


15.7?™ (9.2+|; 8 ) 


C5 


109.61 


97 O+50.1 
z/ - -11.5 


B7 


111.93 


11.5^(6.9^) 


C8 


115.37 


21 6 +26 - 2 

zl - -8.0 








C3 


122.47 


9+10.6 
-5.6 








CD6 


46.26 


5 1+ 13 ' 8 

-3.3 


A2 


24.08 


13 9+5-7 

-3.5 


CD1 andCD2 


70.61 


n 7+17.4 
'■'-4.9 


Al 


96.95 


26 4+ 73 - 6 
zu - -12.2 


CD5 


96.39 


10.2+f., 8 


A3 


108.52 


26 l +58 -° 
ZO -11.2 


CD4 


102.67 


10.5!f/ 








CD3 


142.19 


20.2*5-2 









Table 6: Thermal inertia of an individual scan, r ind , for the base-line cases, derived with fixed 
A v ,o and /f ast ,o- The scans are on the order of the mean value of the phase angle, a. The upper 
limit of Tind is set to be 100 Jm _2 K _1 s 1/2 due to an insufficient coverage in simulations. The 
different combinations of A v ,o and /f astj0 for the high \B'\ case (B1-B4) and the low \B'\ case 
(B5-B8) are used. For comparison, the values of r ind estimated using a single combination of 
A v ,o and / fastj o for all the B ring data are shown in parenthesis. 



Ring 


'"min 


r . max 


1 




(cm) 


(m) 




Cring 


1 


10 


3.1 


B ring 


30 


20 


2.75 


Cassini division 


0.1 


20 


2.75 


A ring (inner) 


30 


20 


2.75 



Table 7: The lower and upper cut-off sizes and the power-law index of ring particles estimated 
in French and Nicholson (2000). The inner A ring is the inside the Encke division. 
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Figure captions 

Fig. 1. Azimuthal temperature variations for the C ring. The title of each panel is the scan 
number, which is on the order of the scan date, whereas the panels are placed on the order 
of the phase angle, a. Green crosses are observed temperatures. Curves are modeled 
temperatures with best- fit parameters for four different cases: (1) a single value of the 
thermal inertia (r ) is used for all particles in all scans (black solid curves), (2) different 
scans can have different values of the thermal inertia r ind while A v and / fast are fixed 
for all scans (blue dashed curves), (3) the thermal inertias for slow and fast rotators have 
different values, r s i ow and r fast , whose values are shared in all scans (red solid curves), and 
(4) the same as (3) except the size of fast rotators, r tast , is 5 mm (purple dotted curves). 
See Table 5 for other fitted parameters. The values of T, r slow , and r fast shown in panels 
are in units of Jm _2 K _1 s~ 1/2 . The values for the reduced chi-square, x 1 , for each scan are 
also shown in parentheses. The mean solar phase angle, a, is shown on the top left of 
each panel. 

Fig. 2. Same as Fig. 1 but for the B ring. The panels are placed on the order of the phase angle 
but separately for the high B' scans (B1-B4) and the low B' scans (B5-B8). The purple 
dotted curve is the case of r fast = 5 mm for the high B' scans and r fast = 1 mm for the low 
B' scans. 

Fig. 3. Azimuthal variations of the temperature (top) and the optical depth averaged over a 
footprint (bottom) for the Cassini division. The purple dotted curve is the case of r fast = 1 
mm. 

Fig. 4. Same as Fig. 1, but for the A ring. The purple dotted curve is the case of r fast = 5 mm. 

Fig. 5. Bolometric Bond albedo, fraction of fast rotators, and thermal inertia for the Saturn's 
main rings for the base-line cases. The same thermal inertia value for all particles is 
assumed in the derivation of A v ,o, /f as t,o> and r (black crosses), whereas different thermal 
inertias for slow and fast rotators are assumed in the derivation of A v (red diamonds), / fast 
(red diamonds), r slow (blue diamonds), and r fast (red triangles). See also Table 5. 

Fig. 6. Contour of the reduced chi-square, &x% sliced on the A v vs. T plane with the best-fit 
value of the fraction of fast rotators, ff ast ,o (from Table 5). Note that the error bars are 
estimated using the surface of A^l = 1.0 in the three dimensional parameter space and 
displaying these sliced contours is only for the purpose of giving some insights on error 
bars at other confidence levels (the surfaces of A^ 2 = 1 .00, 2.7 1 , and 6.63 give confidence 
levels of 68.3 %, 90%, and 99%). 

Fig. 7. Thermal inertia r ind vs. phase angle a for the base-line cases. In the parameter fits, 
Ay (= A v ,o) and /f ast (= /f as t,o) are fixed for all scans in each ring, but the thermal inertia is 
treated as a free parameter. The solid lines represent linear fits of Fi nd from all scans. The 
values of r ind for the C and B rings at a = 4.7 - 6.2° are from Ferrari et al. (2005) (the 
offset in a between these two data is for avoiding overlaps of the marks but ground-based 
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observations for both rings were done simultaneously in two different campaigns). The 
values of r ind are also shown in Table 6. 

Fig. 8. Contour of ring physical temperature on the plane of the vertical and azimuthal lo- 
cations. The vertical coordinate z normalized by the scale height h is positive on the lit 
face. Three cases with different combinations of r and type of vertical motion are shown. 
The used parameters are as follows: the heliocentric distance is 9.1 AU, B' = 21°, r v = 
105,000 km, A v = 0.6, / fast = 1.0, and T = 13 Jm- 2 K- l s~ 1/2 . The vertical dashed lines 
represent shadow lines. 

Fig. 9. Contour of the normalized reduced chi-square, Ax 2 Y , sliced on the r s ] 0W vs. Ff ast plane 
with the best-fit values of A v and /f ast (from Table 5). Note that the error bars are estimated 
in the four dimensional parameter space and displaying these sliced contours is only for 
the purpose of giving some insights on error bars at various confidence levels. 

Fig. 10. Lower limit of size of fast rotators rf ast m i n vs. lower limit of thermal inertia of fast 
rotators r fast min . The size r fast is fixed to be 10 cm for the estimation of r fast min . The 
lines represent theoretical predictions (Eq. Q) for three different illumination cycles, 
In I co. The values of r fast min for the B ring at low B' and the Cassini division are not well 
estimated and assumed to be 1 Jm~ 2 K~ 1 s" 1,/2 . 

Fig. 11. Comparison of fitted parameters between this work and Paper II. The solid curves 
are from Paper II (dotted curves are error sizes) and diamonds are from this work. The 
scale hight ratio h r is assumed to be 3 except h r = 1 for the B ring. The thermal inertia 
of all particles is assumed to be the same. The red lines and marks are for the case with 
standard sinusoidal vertical motion without bouncing whereas the black lines and marks 
are for the case with cycloidal motion which assumes particle bouncing at the midplane. 
The values for the B ring from this work are the mean values of high B! and low B' . The 
radial offset between the red and black diamonds for the B ring is for avoiding overlaps 
of the marks. 

Fig. 12. Thermal inertia as a function of fraction of bouncing particles f\, for the mid B ring. 
The solid lines are the thermal inertia estimated from this study, r , and the dashed lines 
are the thermal inertia estimated in Paper II, r Pn . Two cases of optical depth profiles, Tpps 
and Tuvis> are considered in both works. The value of r PII for the case with r UV is and 
/b = 1 is uncertain due to the limited coverage in parameter fits in Paper II, but at least 
larger than 40 Jm ^K^'s™^ 2 . The offsets in /b are for avoiding overlaps of the marks: the 
simulations were made for either f b = or 1. 

Fig. 13. Thermal inertia vs. fraction of fast rotators for the Saturn's main rings. Black bars 
represent the range of these parameters directly estimated in model fits to observed data 
(To and /f ast ,o)- The red and blue curves are relations between these two parameters pre- 
dicted from Eqs. ( [10] ) - ( 13 ) with the size distribution from French and Nicholson (2000) 



(Table 7). 
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Fig. 1 . Morishima et al. 
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Fig. 1 -continued. 
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Fig. 2. Morishima et al. 
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Fig. 2-continued. 
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Fig. 3. Morishima et al. 
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Fig. 4. Morishima et al. 
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Fig. 6. Morishima et al. 
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Fig. 8. Morishima et al. 
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Thermal inertia vs. fraction of bouncing particles 




0.0 0.2 0.4 0.6 0.8 1.0 

Fraction of bouncing particles 



12. Morishima et al. 



-47- 




13. Morishima et al. 



